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We propose a scheme for preparing an EPR state in posi- 
tion and momentum of a pair of distantly-separated trapped 
atoms. The scheme utilizes the entangled light fields output 
from a nondegenerate optical parametric amplifier. Quantum 
state exchange between these fields and the motional states of 
the trapped atoms is accomplished via interactions in cavity 
QED. 
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I. INTRODUCTION 

In 1935, Einstein, Podolsky and Rosen (EPR) pro- 
posed a now famous gedanken experiment involving a 
system of two particles spatially separated but correlated 
in position and momentum as described by the Wigner 
function 



W(q 1 ,p 1 ;q 2 ,P2) oc <%i + q 2 )S(p 1 -p 2 ), 



(1) 



where q\ and q 2 are the continuous position variables of 
the particles with corresponding conjugate momenta pi 
and pi . With the assumption of local realism, but with 
the apparent ability to assign definite values to canoni- 
cally conjugate variables of one particle from measure- 
ments of the other particle in this system, a conflict with 
the Heisenberg uncertainty principle seemingly follows, 
which led EPR to conclude that quantum mechanics is 
incomplete. Bohm 0] adapted this argument to a system 
of discrete (dichotomic) variables, to which Bell applied 
his classic analysis, deriving the so-called Bell inequalities 
|j|-|5| which quantify explicitly the conflict between lo- 
cal realism and quantum mechanics. Note that although 
measurements of (qi,Pi) do not lead to a violation of a 
Bell inequality for the original EPR state of Eq. (|l|) (Ref. 
j|, p. 196), the entanglement of this state guarantees that 
an appropriate set of variables exists for which a contra- 
diction with local realism would be manifest (|J^] . 

Experimental demonstrations of the conflict between 
quantum mechanics and local realism have concentrated 
almost exclusively on systems of discrete variables, such 
as electron spin or photon polarization HJg]. Only one 
experiment, by Ou et al. following suggestions by 

Reid and Drummond pl| , has in fact realized the EPR 
paradox as originally envisioned by EPR; that is, for 



canonically conjugate variables with a continuous spec- 
trum. The EPR "source" in this experiment was a non- 
degenerate optical parametric amplifier (NOPA) , and the 
relevant variables were the quadrature amplitudes of the 
entangled electromagnetic fields generated in the para- 
metric process. These amplitudes are analogous to the 
position and momentum of a particle and can be mea- 
sured very efficiently via homodyne detection [ |l2| . 

In the present work, we describe a new scheme that 
goes beyond an analogy and actually realizes an EPR 
state in position and momentum for a pair of massive 
particles at distinct physical locations. Our proposal for 
achieving stored entanglement for continuous quantum 
variables is based upon a set of interactions in cavity 
quantum electrodynamics (QED) that allows for the ex- 
change of quantum states between the motion of trapped 
atoms and propagating light fields [jl3|. By exploiting 
these interactions and the light source of [p|Jio| , we show 
that it should be possible to prepare deterministically a 
state of the form ([y) for a pair of trapped atoms located 
at macroscopically-separated sites. Beyond conventional 
(q,p) projections as in homodyne or heterodyne mea- 
surements, the setting of atom traps and cavity QED 
also enables detection strategies for the explicit demon- 
stration of the nonlocal character of the resulting EPR 
state. Moreover, the techniques that we describe could be 
important resources for the realization of quantum net- 
works, a particular example being the creation of EPR 
states to enable the teleportation of the center-of-mass 
wave function of a massive particle |14|. 



II. 



TRAPPED ATOM COUPLED TO AN 
OPTICAL CAVITY MODE 



We begin with the basic setup that facilitates the 
motion-light coupling fundamental to our scheme [13 
this setup was originally considered by Zeng and Lin [15 
We consider a single two-level atom (or ion) confined in 
a harmonic trap located inside an optical cavity. The 
atomic transition of frequency Lu a is coupled to a single 
mode of the cavity field of frequency lu c and is also as- 
sumed to be driven by an external (classical) laser field of 
frequency ■ The physical setup and excitation scheme 
are depicted in Fig. 1. The cavity is aligned along the 
x-axis, while the laser field is incident from a direction in 
the y-z plane (i.e., perpendicular to the a;-axis). 
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The Hamiltonian describing the internal and exter- 
nal atomic degrees of freedom plus the atom-cavity and 
atom-laser couplings takes the form (in a frame rotating 
at the laser frequency for the internal-atomic and cavity 
operators) 



H 



y hvji&tbj + 1/2) + hSdta + hAa + &- 

j=x,y,z 

+ ft [£l(v, z, t)a+ + ££{y, z, t)aJ\ 
+ %o Bm.(kx)(a' <r_ + a+a). 



(2) 



Here, {v x ,Uy,Vz} are the harmonic oscillation frequen- 
cies along the principal axes of the trap, bj and a are 
annihilation operators for the quantized atomic motion 
and cavity field, respectively, <r_ — \g){e\ is the atomic 
lowering operator, and S — uj c — wl and A — w a — lu^. 
The quantity £h(y, z,t) is the (possibly time-dependent) 
amplitude of the laser field. The single-photon atom- 
cavity dipole coupling strength is given by go, while the 
sine function describes the standing wave structure of the 
cavity field (we assume that the centre of the trap is lo- 
cated at a node of the cavity field), with k = 2tt/X the 
wavenumber of the field and x = [h/(2mv x )] 1 ^ 2 (b x + b\,). 

In Jl3|] a number of assumptions and approximations 
are made in order to simplify the model. In particular: 

• The detunings of the light fields from the atomic 
transition frequency are assumed to be very large 
(i.e., A ^> |£l|j<?o> S, Vj), enabling atomic sponta- 
neous emission to be neglected and the internal 
atomic dynamics to be adiabatically eliminated. 

• Any forms of motional decoherence associated with 
the trap itself are ignored. 

• The size of the harmonic trap is assumed to 
be small compared to the optical wavelength 
(Lamb-Dicke regime), enabling the approximations 
sin(fci) ~ rj x (b x -\-b x ), where r\ x (^C 1) is the Lamb- 
Dicke parameter, and £l(j/, z, t) ~ £i J (t)e~' l '^ Ij 



• The cavity and laser fields are tuned so that 5 = 

LU C - UJ L = V x . 

• The trap frequency v x and cavity field decay rate K a 
are assumed to satisfy v x ^> n a » \(g rj x /A)£i,(t)\. 
The first inequality allows a rotating-wave approxi- 
mation to be made with respect to the trap oscilla- 
tion frequency, while the second inequality enables 
an adiabatic elimination of the cavity field mode. 

Given these conditions one can show that the motional 
mode dynamics in the x direction can be described by 
the simple quantum Langevin equation 



b x ~-Y{t)b x + y/2Y{t) 5 in (i) , 



(3) 



where b x = e lVlt b x , Y{t) = [g n x £ L (t) / A] 2 / n a , ando in (£), 
which satisfies the commutation relation [a; n (t), al n (t')] = 



6(t — t'), is the quantum noise operator describing the in- 
put to the cavity field (in a frame rotating at the cavity 
frequency). In this way, the statistics of the input light 
field can be "written onto" the state of the oscillator. 
In E| it was shown how this effect can be used to ef- 
ficiently prepare a squeezed state of the motion of the 
trapped atom. Here we extend that work further to the 
generation of entanglement between the motional states 
of trapped atoms at separated sites. In particular, our 
protocol transfers entanglement from a pair of quantum- 
correlated light fields to a pair of trapped atoms in a 
process of quantum state exchange, or qusex. 



III. LIGHT SOURCE: NONDEGENERATE 
PARAMETRIC AMPLIFIER 

Our source of quantum-correlated light fields is taken 
to be a NOPA operating below threshold The 
light fields may be nondegenerate in polarization or in 
frequency. We denote the annihilation operators for the 
two intracavity field modes, of frequencies u)\ and ci>2> by 
ci and C2, and the Hamiltonian describing the coupling 
between these modes takes the form (in a rotating frame) 



Hi = ihe [c\C2 



at 



(4) 



where e is the coupling strength, proportional to the non- 
linear susceptibility of the intracavity medium and to 
the strength of the coherent pump field (at frequency 
OJl + lj 2 ). 

Assuming the cavity mode amplitudes to be damped 
at the same rate k c , equations of motion for the mode 
operators (in the rotating frame) can be derived as 



Cl,2 



( 1 2) 

where c in ' it) are the vacuum input fields to the NOPA 
cavity modes (see, e.g., |l6| , |l7f ). The output fields from 
the NOPA then follow from the boundary conditions 



-k c ci,2 - ec 2 1 



s(i.2) 



(t), 



(5) 



4t 2) w+4:' Zj w = ^c 1:2 w, 



s(1.2) 



(6) 



The (linear) equations above are readily solved in 
a Fourier-transformed space defined by Z{uj) = 
(27T)- 1 / 2 / dtZ(t)e~ tuJt . Defi ning quadrature phase am- 
plitudes ("positions" and "momenta") for the output 
fields by 



X, 



(1.2) 
out 

(1,2) 



(t) 



a(l,2) 



s(l,2)t 



(t), 



(t) = -i{&(t)-&\t)} , 
the sum of the X amplitudes is derived as |H| 



(J) 
(8) 



out V 

K r -e 



K c + e — ILO 

as 6 — ► k 



and u> 



X 



(2) 



o, 



(9) 
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while the difference of the Y amplitudes is 



This master equation has a steady state solution 



y«( w )-yp)( w ) 



IUJ 



k c + e — iuj 
as e — > k c and u> 



0. (10) 



So, the two output fields are highly correlated and, close 
to uj = and for e — > k c , their quadrature amplitudes 
exhibit precisely the properties of the original EPR state, 
as demonstrated explicit ly by the Wigner function for the 
state of the (1,2) fields 



t ry by t 

00 



IV. LIGHT-TO-MOTION QUANTUM STATE 
EXCHANGE 



As depicted in Fig. 2, the two NOPA output fields are 
assumed to be incident on separate cavities, each contain- 
ing a trapped atom in the configuration described earlier. 
Note that the output fields from the NOPA are resonant 
with the respective cavity mode frequencies. We assume 
that r(i) — r, a constant, and, for simplicity, that T is 
the same for both configurations. Denoting the motional 
mode operators for the two atoms along the x-axis by b\ x 
and b 2x , respectively, the two systems are thus described 
by fliT 



bj x J^bj x 
= -Tb jx 



5« 



2TQt-^ (i = 1, 2) 



(11) 



where r is a time delay (assumed the same for both cav- 
ities); provided the coupling between the NOPA and the 
cavities is unidirectional, this delay can essentially be ig- 
nored OUTT 



If the bandwidths of the input light fields from the 
NOPA are sufficiently broad, in particular if k c 3> T 
[i.e., c^t\t) can be regarded as quantum white noise 
operators in (jll|)1, then one can perform an average over 
the input fields and derive a master equation (see, e.g., 
Jl6],[l7|) for the density operator p of the motional modes 
alone, 



p= T(N + l)(2b lx pb\ x - b\ x b lx p - pb\ x b lx ) 
+ TN(2b\ xP b lx - b lx b\ xP - pb lx b\ x ) 
+ T(N + l)(2b 2x pP 2x - b\ x b 2x p - pb\ x b 2x ) 
+ TN(2bl x pb 2x - b 2x b\ xP - pb 2x b\ x ) 
+ 2TM(bi x pb 2x + b 2x pbi x - bi x b 2x p ■ 



- 2rM (b\ x pbl + b\ x pb\ x - b\ x b\ x p - pb\ x b\ x ) , (12) 

with the parameters iV and M given, in terms of the 
NOPA parameters, by 



pb\ x b 2x ) 



N = 



2„2 



M = 2k c c 



(Kl e 2 ) 2 



(13) 



P ss = IViaX^al, (14) 
i.e., a pure state, with 

l^u) =Sia(r)|0)i B |0)2x 

OO 

= [cosh(r)] -1 ^2 [-tanh(r)]" 1 \m) lx \m) 2x , (15) 

m— 

where \m)i Xj2x are Fock states of the motional modes and 
S 12(f) is the two- mode squeezing operator 



Si 2 (r) = exp 



r [bi x b 2x 



b\J>L 



(16) 



This operator transforms the mode annihilation opera- 
tors for the atomic motion as 



S\ 2 (r)bi x Si 2 (r) = cosh(r) b lx - sinh(r) b\ x , 
S\ 2 (r)b 2x Si 2 (r) = cosh(r) b 2x - sinh(r) b\ x , 



where cosh(r) = \/N + 1 and sinh(r) 
position and momentum operators as 



it follows that 



P, 



-i b 



S\ 2 (f) (Qi + Q 2 ) Si 2 (r) = e~ r (Qi + Q 2 ) 
S\ 2 (f) (Pi - P2) Si 2 (r) = e~ r (Pi - P 2 ) , 



(17) 
(18) 

N. Defining 
(19) 



(20) 
(21) 



and so, in the limit e — > k c (i.e., r 00), an EPR state 
in the positions and momenta of the two trapped atoms 
is established. 

The nature of the correlations inherent in the joint 
state ( |l5|) of the atomic motion is most clearly expressed 
through the Wigner function for this state : 



W(qi,pr,q 2 ,p 2 ) 
4 
72 



in)' e'-'} 



exp{- [(qi +q 2 f + (pi 
X c^{-[(qi-q 2 ) 2 + (pi+p 2 ) 2 }e- 2r } (22) 



2 I „+2r1 

2 ] e-*r) 

C S(qi + q 2 )S(pi -p 2 ) as r ^ 00 , (23) 



with C a constant. This entangled state is achieved in 
steady state over a time t ^> The coupling to the 

external fields from the NOPA can then be turned off by 
setting £li,2 to zero. The result is a stored EPR state for 
the motion of two trapped atoms that would persist for a 
duration set by the timescale for motional decoherence. 



V. DISCUSSION 

Before considering some of the interesting possibilities 
offered by this system, we return briefly to some of the 
major assumptions associated with the model. Firstly, 
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the finite effect of atomic spontaneous emission events 
on the motional state can be estimated as in 13 1. This 
effect can be neglected provided the rate of these events 
is much smaller than the rate T at which the motional 
steady state is achieved. The condition one derives by en- 
forcing this inequality amounts to the condition of a "one- 
dimensional" atom in cavity QED, C\ = <7o/(fta7) ^S> 1, 
where 7 is the linewidth (FWHM) of the atomic tran- 
sition |20[ ]. Note that the experiment of Ref. j2lj has 
achieved C\ — 70. With regards to the trapping po- 
tential, harmonic frequencies on the order of tens of 
MHz have been achieved in ion traps, with corresponding 
Lamb-Dickc parameters on the order of 0.1 and smaller 
p3| . Note, however, that large values of the entangle- 
ment parameter r imply population of large-m number 
states and a broader spread of the atomic wavepacket. 
Given that the mean excitation number for the state (|l5| ) 
is n = sinh 2 (r), a more precise form of the Lamb-Dicke 
assumption would be rj x \Jn + 1 = rj x cosh(r) -C 1. With 
trap frequencies such as those quoted above, the condi- 
tion v x 3> K a should be satisfied for a cavity field decay 
rate of a few MHz or less; such values of K a have not 
been realized in current experiments as in Ref. pl| , but 
potentially could be with improved cavity finesse as in 
Ref. p^j . Assuming that this is the case, likely magni- 
tudes for the rate T would then be tens or hundreds of 
kHz. Finally, timescales for motional decoherence and 
heating in recent ion trap experiments are of the order of 
milliseconds, with further improvement likely p3[ ; given 
the various rates discussed above, these effects would not 
be expected to hamper the preparation of the entangled 
state. 

As for applications of this system, further investiga- 
tion of the EPR paradox would obviously be possible, 
with a variety of motional state measurements able to 
be implemented on the trapped atoms [ p4|j25| , possibly 
also via the cavity field In particular, violations of 
a Bell inequality for the state (051) can be obtained with 
measurements that project onto a basis of even and odd 
phonon number for each of the trapped atoms ^,0] . To 
the extent that a "macroscopic" number of quanta may 
in principle be involved, such investigation could also ad- 
dress new viewpoints on the compatibility of quantum 
mechanics with local realism p6| . 

On a somewhat more applied side is the possibility of 
using the EPR state ([l5]) for quantum dense coding p7| , 
or for the teleportation of the quantum state of a system 
with continuous variables [Bq , B9[ , generalizing the origi- 
nal discrete- variable teleportation protocol of Bennett et 
al. J|(| . This elegant adaptation of the EPR paradox has 
in fact been realized with light fields, again using opti- 
cal parametric amplifiers and homodyne measurements of 
quadrature amplitudes | pl| . The scheme outlined in this 
paper opens the door to teleporting an atomic center- 
of-mass wave function fl4| , by providing the motional 
state entanglement required by the continuous-variable 
teleportation protocol. 

Such a capability is also of considerable interest in the 



related context of quantum computation with trapped 
atoms and light f32| , |33|| . Here, we specifically have in 
mind protocols that combine quantum information pro- 
cessing with both discrete and continuous variables [p4| . 
Any implementation of a qubit (e.g., internal atomic 
states or photon polarization) could be linked with an 
external degree of freedom (e.g., atomic center-of-mass 
or complex amplitude of the electromagnetic field) , with 
the complete system viewed as a composite unit (qubit 
plus qunat |34|]) for protocols such as quantum commu- 
nication between distant nodes of a quantum network 

mm- 

Further to this theme, note also that we need not re- 
strict ourselves to a single trapped atom at each site. For 
example, if there are K atoms inside each cavity, then fo- 
cussing the coupling lasers (£1,1,2) sequentially on atoms 
1, 2, . . ., K at each site (and neglecting any direct inter- 
action between neighbouring atoms at each site) would 
generate a set of pair-wise EPR-entangled atoms. Al- 
ternatively, and perhaps more interestingly, one might 
consider the case in which the K atoms at each site are 
simultaneously coupled to the cavity. Assuming for sim- 
plicity that they have identical coupling strengths, then 
the system would again be described by Eq. (|jlj|), but 



iU) 



1 u l,2x 



with the replacements 61.2a; — ► -81,22; = K^ 1 ' 2 
and r — > KT. Such dynamics evidently leads to a highly 
entangled state of all 2K atoms, a situation of potentially 
great utility and, indeed, of considerable general interest. 
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FIG. 1. Schematic of proposed experimental setup and ex- 
citation scheme for coupling between the motion of a trapped 
atom and a quantized cavity mode of the electromagnetic 
field, and thence to a freely propagating external field. Note 
that all input and output to the atom-cavity system is as- 
sumed to be through just one mirror; the other mirror is as- 
sumed to be perfect. 



FIG. 2. Preparation of an EPR state of the motion of 
two separated trapped atoms. The output modes from a non- 
degenerate parametric amplifier (NOPA) are incident on two 
separated trapped-atom-cavity configurations of the type de- 
picted in Fig. 1. Faraday isolators (F) facilitate a unidirec- 
tional coupling between the entangled light source and the 
atom-cavity systems. Note that in practice the outputs from 
the two NOPA modes might actually be through the same 
mirror, but could be separated due to differing polarizations 
or frequencies. 
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